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We investigate the conductance through and the spectrum of ballistic chaotic quantum dots 
attached to two s-wave superconductors, as a function of the phase difference cj> between the two 
order parameters. A combination of analytical techniques - random matrix theory, Nazarov's circuit 
theory and the trajectory-based semiclassical theory - allows us to explore the quantum-to-classical 
crossover in detail. When the superconductors are not phase-biased, cj> = 0, we recover known results 
that the spectrum of the quantum dot exhibits an excitation gap, while the conductance across two 
normal leads carrying JVn channels and connected to the dot via tunnel contacts of transparency 
Tn is cx Fn-ZVn. In contrast, when (j> = n, the excitation gap closes and the conductance becomes 
G oc TnA^n in the universal regime. For Tn < 1, we observe an order-of-magnitude enhancement 
of the conductance towards G oc Nh in the short-wavelength limit. We relate this enhancement 
to resonant tunneling through a macroscopic number of levels close to the Fermi energy. Our 
predictions are corroborated by numerical simulations. 

PACS numbers: 74.45.+C, 73.23.-b, 74.78.Na, 05.45. Mt 



I. INTRODUCTION 

Andrccv reflection plays a central role in the descrip- 
tion of hybrid nanostructures with both superconducting 
and normal components [l[ . The process is illustrated in 
Fig. [TJ When a negatively charged electron in the nor- 
mal metal hits the interface with a superconductor, it is 
retroreflected into a positively charged hole. The super- 
conductor acts as a charge-inverting mirror, which leads 
to a variety of intriguing effects in transport and spec- 
troscopy of hybrid nanostructures 0. Here, we show 
how Andreev reflection resonantly enhances the tunnel- 
ing conductance through a quantum dot by orders of 
magnitude, when the dot is coupled to two supercon- 
ductors with phase difference <fi = n between their order 
parameters. The effect quickly disappears as (j) is tuned 
away from resonance and might have applications in low- 
temperature current switching devices and magnetic flux 
"transistors" , where (f> would be controlled by an external 
magnetic flux. 

In a chaotic Andrccv billiard, an impurity-free quan- 
tum dot coupled to one or more superconductors, An- 
dreev reflection renders all orbits periodic near the Fermi 
level Q . Periodic orbits consist of an electron and a hole 
retracing each other's path, and therefore touching the 
superconductor at both ends. This fact suggests that 
the density of states of an Andreev billiard can be de- 
termined by Bohr-Sommerfcld quantization. For a single 
superconducting lead, this results in a density of states 
which is exponentially suppressed for low energies, on the 
scale of the Thouless energy Et,s = ^/2td,s, where td,s 
is the average time between Andreev reflections [1, Q . Si- 
multaneously, one expects random matrix theory (RMT) 
to be valid for chaotic quantum dots RMT however 
leads to a different prediction, that a hard gap in the 




Figure 1: Andreev reflection by a superconductor (S). An 
incoming electron (e) in a normal metal (N) near the Fermi 
energy Ef is reflected as a hole (h) with opposite velocity. 

density of states opens up at an energy Eq ~ 0.6£t,s [li- 
lt was realized in Ref. [7j that the discrepancy between 
the theories was not a short-coming in one of them but 
indicates that they have different regimes of validity. The 
crossover between the two regimes is determined by the 
Ehrcnfest time, which is the time scale for an initially 
minimal wave packet to spread over the system. For a 
closed chaotic billiard of linear size L c and Lyapunov ex- 
ponent A, the Ehrcnfest time reads te = A^ 1 ln[/cpic], 
with the Fermi wavenumber kp of the quasiparticles Q . 

RMT is valid in the universal regime of vanishing 
Ehrcnfest time, te <C td,s> while Bohr-Sommcrfeld quan- 
tization gives a good description in the opposite semi- 
classical regime, te 2> ttj.s- Previous investigations of 
the crossover from the RMT regime to the semiclassi- 
cal regime in Andreev billiards have considered the den- 
sity of states in the presence of a single s-wave super- 
conductor, focusing on the reduction of the excitation 
gap in the semiclassical regime, as well as its fluctua- 
tions @, i, m m m, m ei m, m, 53, m m. The 

Ehrcnfest time effects arc usually small and it takes a 
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large computational effort to extract them numerically. 

In this article we investigate transport via normal 
metallic leads through an Andreev billiard connected to 
two 0-biased superconductors. The system we consider 
is shown in Fig. [5] A chaotic quantum dot with mean 
level spacing S is coupled to two superconductors. The 
superconducting lead 1 (2) has order parameter Aoe 1 * 1 
(Aoe 1 ^ 2 ), and we take Ao £ 3?. One way to control the 
phase difference (f> = 4>\ — 02 is to have them form a 
loop, which one then threads with a magnetic flux 4> - 
this is sketched in Fig. [2] Alternatively, the phase dif- 
ference can be generated by a supcrcurrent. In this case 
<p = 4>\ — 4>2 = 27r$/<l>0; i n terms of the flux quantum 
$0 = h/2e. Since any global phase can be gauged out, 
we set (j>i — cf>/2 and 02 = —0/2. We assume that each 
superconducting lead has Ns channels connected to the 
quantum dot via tunnel contacts with tunnel probability 
Tg. The average time td,s = ^/2-Et,S between Andreev 
reflections is associated with the superconducting Thou- 
less energy i?T,s = 2Asrs^/47r. Quasiparticle excitations 
have energy E measured from the Fermi energy E-p , with 
\E\ <C Ao. Thus the quasiparticlcs cannot penetrate into 
the superconductors. We also assume Ao <§C Ep, in which 
case Andreev reflection perfectly retroreflects electrons 
into holes and vice versa, with a phase shift — 7r/2 ± 0j. 
At = 0, this Andreev phase shift is solely due to the 
penetration of the wavefunction into the superconduc- 
tor. The additional phase shift ±0i is due to the global 
phase of the superconductor at contact i, where the mi- 
nus sign is for reflection from electron to hole, and the 
plus sign for reflection from hole to electron. We will re- 
fer to the quantum dot with two superconducting leads 
(but no normal lead) as the closed Andreev billiard. For 
the transport set-up we connect the quantum dot to two 
external normal leads through a tunnel contact with tun- 
nel probability Tn- The current is measured between the 
right (R) and left (L) lead, which both carry An trans- 
port channels. 

The combination of Andreev and Ehrenfest physics in 
such an Andreev interferometer gives rise to an ordcr- 
of-magnitude enhancement of the tunneling conductance 
when = 7r. In the tunneling regime Tn — > 0, the con- 
ductance is strongly affected by the density of states of 
the closed Andreev billiard. In the absence of supercon- 
ductivity (Ts = 0) and for TnAn > 1 (we restrict our- 
selves to that regime to avoid Coulomb blockade effects) , 
the conductance is just the classical series conductance 
G = G r N 7V N /2 + 0(1), where G = 2e 2 /h $2%. In 
contrast, we will see below that in the limit te/td ~ 1, 
and for well coupled superconductors, Tg — 1, the con- 
ductance at = tt becomes G oc Go An- We iden- 
tify the mechanism behind this enhancement as macro- 
scopic resonant tunneling through a large number of 
low-energy quasi-degenerate Andreev levels very close to 
the Fermi energy. Indeed, for <f) = tt, all periodic tra- 
jectories touching both superconductors, irrespective of 
their lengths, contribute to the density of states at the 
Fermi energy in the semiclassical regime. Such a large 




Figure 2: Quantum dot doubly coupled to a superconducting 
ring. Both contact to the superconductor have Ns channels 
and tunneling probability Ts. The flux threading the ring 
causes a phase difference (j> = 27r$/"l>o between the two su- 
perconducting contacts. Transport through the structure is 
investigated between two external normal leads (L and R), 
each carrying JVn transport channels, with tunnel contact of 
tunnel probability Tn- 



effect is already clearly visible for moderate Ehrenfest 
time, in contrast to previous works which found rather 
small manifestations of Ehrenfest physics for larger te/td 
03 EL OH El H, [H, El. While similar behaviors 
were reported for cavities without transport (2l| or inter- 
nal [22| mode-mixing, this effect was not found in earlier 
analytical works on transport through chaotic Andreev 
interferometers 0, |M HE HI • 

The paper is organized as follows. In section II we cal- 
culate the density of states of the closed Andreev billiard. 
We repeat previous results in the RMT limit and calcu- 
late the Bohr-Sommcrfcld density of states as a function 
of 0. We show that due to a special resonance condition, 
the density of states for phase differences close to tt has a 
large peak around the Fermi energy. We proceed with a 
calculation of the conductance in Sect ion [Hll In the RMT 
limit we use Nazarov's circuit theory 2?], 28| ■ For weakly 
coupled normal leads, we find that the conductance is 
small, cx T^An, for <fi <C tt, while it reaches the classical 
series conductance for (f> = tt, because of the disappear- 
ance of the gap in the density of states. In the semi- 
classical limit we emplo y a semiclassical trajectory-based 

method [U, M HE OH 1H ■ We find that the lar g e P eak 
in the density of states at the Fermi level has its coun- 
terpart in the conductance, which is a factor 1/Tn 3> 1 
larger than in the RMT limit. In section IV, numerical 
simulations confirm the validity of our analytical results. 



II. SPECTRUM OF ANDREEV BILLIARDS 

A. Random matrix theory 

The density of states in the universal regime has been 
calculated in Refs. [1, [13 using RMT. We repeat the 
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main steps of the calculation for Ts = 1, and refer the 
reader to Refs. [1, [34| for the calculation at arbitrary 
rg. At low excitation energy £ « Ao, the spectrum of 
Andreev billiards is most conveniently obtained from the 
low-energy effective Hamiltonian [34[ 

or _ ( H -nWW T \ n s 

* ~ \ -ttWW t -H* ) ' 

rrr _ , (MS\^ 

\ir 2 J 

m=l,2, ...M, 7i = l,2,. 



with Bcsscl functions Jo and J\. This prediction retains 
its validity down to low excitations E < 6. For energies 
E > 6 (but still E < A ), Eq. J5]) reduces to the mean 
field theory and gives p(E) = 2/5, while for E < 5 it 
predicts a linearly vanishing density of states p(E) = 

TT 2 E/6. 



■2N S , 



(lb) 

B. Semiclassical regime 



which is constructed from the quantum dot's onc- 
quasiparticle M x M Hamiltonian matrix H and the 
projection matrix WW T giving the coupling between 
the quantum dot and the superconductors. We define 
a 2 x 2 Green function G(z) from the matrix Green func- 
tion g(z) = ((z- ■K)- 1 ) 

G= (Gu G iA=i( Tr Sn TrS 12 \ 
\ G 2 i G22 J 7T \ TrS2i TrS22 / ' 

from which the density of states reads 

Prmt(E) = -^Im[G 11 (£; + i + )]. (3) 


Here, G\\ is self-consistently determined by the equations 



The retroreflection at the superconductor makes all 
classical trajectories near the Fermi level periodic [||. 
The elements of a periodic trajectory are an electron and 
a hole segment retracing each other, and touching a su- 
perconducting contact at both ends. The phase accumu- 
lated along a periodic orbit of period T in one period 
consists of (i) the phase 2ET/H acquired during the mo- 
tion in the normal region and (ii) the phase shift acquired 
at the two Andreev reflections. This shift equals — ir if 
the two reflections take place at the same superconduct- 
ing lead and — 7r ± cf> if they take place at different su- 
perconductors. Summing the different contributions and 
requiring that the phase accumulated in one period is a 
multiple of 27r leads to the Bohr-Sommerfeld quantiza- 
tion conditions 
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G U {Z): 

G\ 2 {z) =l + G\ 1 {z) 



G l2 {z) \ 
N S S V 121 ' cos (0/2) ) 



For = 0, the density of states has a gap at an energy 
Eq w 0.6-Et,S- The energy of the gap is reduced when 4> 
is increased, closing at a value <p = it, where one finds a 
flat density of states p(E) = 2/8. 

The theory of Refs. [j, [l7| is only valid for energies 
larger than the level spacing 5. For a phase difference 
<f> = 7r, the system belongs to the symmetry class CI of 
Ref. [35|, and the density of states reads 



( 4a ) E n = Kh(n + 1/2) j T, (6a) 

( 4b ) E n = Trh(n + 1/2 + 4>/2it) j T, (6b) 

E n = nh(n + 1/2 - 0/2tt) / T. (6c) 



P{E) 



2nE/5 



dtJ (t)Jl(t)/t, 



Eq. (|6ap applies to trajectories touching the same super- 
conductor twice, while Eqs. (|6b[) and (|6c|) to trajectories 
touching both superconductors, depending on whether 
the Andreev reflection from electron to hole takes place 
at the superconductor with phase <p/2 [Eq. (|6b[) ] or —0/2 
[Eq. ([6c]) ]. Assuming ergodicity, the corresponding tra- 
jectories have a weight 1/2, 1/4, 1/4, leading to a mean 
W density of states (see also Ref. [36| for a single supercon- 
ducting contact) 
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2tt j T 



-6(E-[r 



TrS' 

2tP T\ 



(7) 

Here P(T) is the classical distribution of return times to a superconducting contact. Inserting the chaotic distribution 
P(T) = cxp (-T/t d ,s)/td,s into Eq. © gives 
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coth 
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(8) 
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with u = E/Et,s- For = Eq. ([8]) reduces to the Bohr-Sommerfeld approximation of Rcfs. 0, H[ and the density 
of states is exponentially suppressed at low energies. 
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Figure 3: Comparison of the density of states calculated 
from Bohr-Sommerfeld quantization and random matrix the- 
ory (RMT). The top panel is for a phase difference <f> = 0, 
and the bottom panel for <f) tt. 



C. Comparison of the two regimes 
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Figure 4: Circuit for our system. We have two normal elec- 
trodes (1,2), two superconducting electrodes (3,4) and one 
internal junction, the chaotic cavity (5). The superconduct- 
ing (normal) leads have Ns (A^n) channels with transmission 
probability Ts (Fn). 



Fig. [3J illustrates the predictions of the previous two 
paragraphs. In the RMT limit we combine the solution of 
Eqs. (|4aj) and (|4bj) with Eq. ((3]), while in the scmiclassical 
limit we use Eq. ((?]). For = (top panel), the density 
of states in the RMT and semiclassical limit are very 
similar. Both densities of states are suppressed at the 
Fermi energy, while they are restored to a value of 2/5 
at an energy scale set by -Et.s- The main difference is 
in the way that p is suppressed: a hard gap in the RMT 
limit vs. an exponential suppression for the semiclassical 
density of states. 

In contrast, when <j> sw tt, the difference between RMT 
and Bohr-Sommerfeld predictions is huge. While in the 
RMT limit the density of states is flat [up to the minigap 
of order S predicted by Eq.([5]) which is not resolved in 
Fig. [5] , it shows a large peak around the Fermi energy in 
the semiclassical limit. For = tt the peak is described 
by a Dirac 5- function (cf. Eq. (flT}|) ). which is why we 
show a density of states at a phase difference slightly 
different from tt in Fig. [31 



The low energy approximation of Eq. ([5]) reads 

2tt e-^-^l u {ix - (j)) 



lim pbs(u) = 

u— >0 



2u 2 



(9) 



For 7T — <p -C 1 this gives a sharply peaked function with 
a maximum value of pss = 
In the limit d> —> tt wc find 



e 2 (ir- 



at energy u = ^-j^- 



lim lim pbs{u) 

U—>0 (b—>TT 



2tt 



-S(u). 



(10) 



Eq. (fTTJ|) predicts a number 0[ATg] of levels around E = 0. 
The peculiar behavior at <fi = tt is due to the fact that the 
two quantization conditions (|6b[) and (|6c[) have solutions 
E = simultaneously for all trajectories, irrespective of 
their period. 



III. TRANSPORT 

A. Random matrix theory 

We use Nazarov's circuit theory of Andreev conduc- 
tance [27| to calculate the two-terminal conductance 
through the system shown in Fig. [5] (see also Ref. [HI). 
The theory describes hybrid nanostructures in the same 
language as electrical circuits: i.e. their structure is re- 
duced to internal junctions, normal/ superconducting ex- 
ternal junctions (electrodes) and leads connecting the dif- 
ferent elements. The leads are characterized by a set of 
transmission eigenvalues. The conductance between dif- 
ferent electrodes is determined by a simple set of alge- 
braic equations. As the energy dependence of the scat- 
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tcring matrix is neglected, the theory is valid for small 
temperatures and voltages. The superconductors have 
the same voltage. The results are valid to leading order 

in r s A s , r N iv N . 

Our system consists of two normal external junctions 
(denoted by indices 1,2), two superconducting external 
junctions 3, 4 and one internal junction, the chaotic cav- 
ity 5. This is illustrated in Fig. 2J The calculation of 
the conductance starts with the assignment of a com- 
plex amplitude to each junction. The amplitudes 
are somewhat analogous to voltages: if a junction n is 
connected to a junction m by a lead with a large con- 
ductance, then /„ — ► f m . The physical meaning of the 
amplitudes is the following: if electrons are injected into 
the k th junction via an additional weakly coupled normal 
electrode, a beam of holes with intensity \fk\ 2 times that 
of the electron beam is retro-reflected into that probe. 
For the external junctions the value of the amplitudes is 
prescribed by the theory: for the normal external junc- 
tions /i,/a = and for the superconducting junctions 
fs = exp {—i(j)/2), fi = exp (i(f>/2). The amplitude f$ can 
be found by the spectral current conservation rule [28| . 
One has /s = tan(#5/2), with #5 to be determined by 



rvvg 
r N 7v N 



cos^ cos6» 5 (l + cos6> 5 ) I 1 + (1 - T N ) tan 2 



= sin0 5 (2-r s +rssin0 5 cosp. (11) 

Once the complex amplitudes are known, the RMT- 
averaged conductance between junctions 1 and 2 is given 
by 



(G)/G 



r N iV N cos6> 5 + T N sin 2 [(9 5 /2] 
2 (l-r N sin 2 [0 5 /2]) 2 ■ 



(12) 



The self-consistency equation pT|) for #5 combined with 
Eq. (fl"2)) determines the conductance. For ballistic leads 
(Tn = Tg = 1) we reproduce the results of Ref. [25j | . 

Of most interest to us is the regime Tn <C 1 and Tg = 
1, where we expect a strong connection between conduc- 
tance and density of states. In the limit 8(f) = ir—tft ^ IV , 
we use an expansion in terms of Tn , while in the opposite 
limit, we first expand in the small parameter 8(f> . From 
Eqs. (fTTj) and (fl"2"|) we find, to leading order in Fn, 



(G)/G 



T 2 N N N ( 1 + 



jy N l+cos ■ 

Ns cos ? 



r N iV N /2 + 0(<5</> 2 



'4, if 50 » r N , 

if 5(f> < r N . 

(13) 



The main feature of Eq. (fT3|) is that, for (f> = it, one 
has the same conductance as a two-terminal quantum 
dot without superconductor. This result agrees with the 
density of states where the influence of the superconduc- 
tors also disappears and the density of states becomes 
flat (bottom panel of Fig. [3]). The minigap of magni- 
tude 8 is hardly resolved, due to the coupling to external 



leads and the associated broadening of levels. On the 
other hand, for (f> ^ tt the density of states is strongly 
suppressed at the Fermi level (top panel of Fig. [3]) and 
we expect a conductance which is reduced due to the 
presence of the superconductors. Eq. (fT3)) does indeed 
predict a conductance of order T 2 ^. 

From Eq. (fT3"|) it follows that the crossover between the 
two limits occurs when 8(f) — T-^N-^/Ns- We now con- 
nect this to the density of states of the Andreev billiard. 
Taylor expanding Eqs. ([3]), (|4"aj) and (|4b"|) in 8(f) gives a 
gap in the density of states at an energy 

iV s 8 

E =8(/)E T , s =S ( f)^. (14) 

Z7T 

The density of states of the closed cavity is broadened due 
to the coupling to normal leads. The coupling to leads re- 
sults in a non-hermitian Hamiltonian H — > H—iirW^W^, 
with H the Hamiltonian of the closed system and Wn an 
M x 2JVn matrix describing the coupling to the normal 
leads H3,H3|. The 2iV N x 2N N matrix W^W N has eigen- 
values w — T^MS/An 2 + O(r^). For the broadening we 
need the eigenvalues of the M x M matrix WnWjJj. As- 
suming that all levels in the dot are coupled to the leads 
in the same manner, we multiply w by 2N^/M to esti- 
mate the broadening as 8E = 2ttwN-^/M = T-^N-^8/2tt. 
Comparing this with Eq. Q14j) we find that 8E = E n for 
Tn^Vn = Ns8(f>. This gives the crossover energy scale for 
Eqs. (fl~3|) . We conclude that the conductance is not af- 
fected by the presence of the superconductors once the 
gap in the broadened density of states closes. 

We finally comment on the effect of a finite tunnel 
barrier with Tg < 1 between the superconductor and 
normal metal. For Tn 1, Eq. (fl"2")) becomes 



l+cos ■ 



'4, if50>r N , 

(G)/G 

«r N . 

(15) 

Here, iVg.eff = cTgA^g where c £ (|, 1) is a numerical fac- 
tor depending only weakly on (f> and Tg . A finite Tg thus 
reduces the number of effective superconducting channels 

to A^S.eff- 



B. Semiclassical regime 

The starting point of our semiclassical treatment is the 
two-terminal conductance through a metallic system in 
contact with superconductors (see Fig. [5]) [38j |. 

G/Go = + + 2 ; tlf ^ ■ (16) 
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We used the transmission probability for a quasi-particle 
of type a = e,h from the normal lead i — L, R to a quasi- 
particle of type (3 = e, h into the normal lead j = L, R, 



T 



0a 



J 1 



T Pa .1 



(17) 
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Here, T^". m i gives the transmission amplitude from 
channel m in lead i to channel n in lead j. Since we 
consider a symmetric configuration, where both normal 
leads carry the same number N-^ of channels, all of them 
coupled to the cavity with the same tunnel probability 
Tn, Eq. (fl"6|) simplifies to 



(G)/Go = (T&) + (Til) 



(18) 



Eq. (|18|) gives the ensemble-averaged conductance to 
leading order in Tn iVjy ^S> 1. 




Figure 5: Schematic of the class I (blue) and class II (red) 
periodic orbits of the closed Andreev billiard, giving rise to 
macroscopic resonant tunneling at (f> = 7r when the cavity is 
connected to external electrodes (L and R). 



To evaluate the resonant contributions to (T^ L ) and 

ssical approach of Rcf 
) Scmiclassically, the trans- 



(T^£), we follow the semiclassical approach of Ref. [2 



L LL 

(see also Refs. [301, L3JJ, [32J, |3 

mission amplitudes read, 



[3a 



= -(2-Kih) 



-1/2 



dj/ ( 



x(^|y)(2/o|m) exp[iS 7 /7i + i7r^ 7 /2] , (19) 

where |m) is the transverse wavefunction of the m th lead 
mode. This expression sums over all trajectories 7 (with 
classical action S 7 and Maslov index /x 7 ) starting at yo 
on a cross-section of the injection (i) lead and ending 
at y on the exit (j) lead, with even (a = (3) or odd 
(a 7^ f3) number of Andreev reflections. The transmis- 
sion probabilities are then given by a double sum over 
trajectories. After the semiclassical approximation that 
£„(y»Hj/> - <%' - V) M one has 



T 



/3a 



dy' £ A yl A; 2 cMiSS/h]. (20) 
71,72 



convert an a quasiparticlc into a (3 quasiparticlc. The 
phase 5S = S 1 \ — S 7 2 gives the difference in action phase 
accumulated along 7I and 72. In the presence of tunnel 
barriers, the stability amplitudes A 1 are given by [HI, [39[ 



_Z> 7 tjtj 



n 



'7v 



M*0 



(21) 



where ly(k) gives the number of times that 7 is re- 
flected back into the system from the tunnel barrier at 
k = L, R, Si, and S2, the transmission and reflection am- 
plitudes satisfy \ti\ 2 = (1 — |r,| 2 ) = r NjS (fori = L, R, Si, 
or S2), and B 2 = (dp yo / dy) 1 measures the rate of change 
of the initial momentum p ya as the exit position y of 7 
is changed, for a fixed sequence of transmissions and re- 
flections at the tunnel barriers. 

There are four types of trajectories to consider: 

(i) paths that do not touch the superconductors, 

(ii) paths that hit one of the superconductor once, 

(iii) paths that hit the same superconductor twice, 

(iv) paths that hit both superconductors. 

Trajectories of type (iv) are flux-dependent. They are 
depicted in Fig. [SJ and it turns out that they are the 
ones giving rise to macroscopic resonant tunneling. Our 
semiclassical investigations focus on those paths. 

We use Eqs. (|2U)) and (|21[) to evaluate the dominant 
semiclassical contributions to Eq. (fT5)) . We subdivide the 
relevant trajectories into class I trajectories, contributing 
to (T^£) (blue trajectory on FigEJ), and class II trajec- 
tories, contributing to (T^ L ) (red trajectory on FigO. 

Class I trajectories are made of the following sequence 



(p) (e) , (h) , (h) . (e) 

ll =Tsl +TsX +P>< Ts3 +Ts3 



(e) 

7si 



(h) 



(22) 



where si and s3 can be interchanged, and p = 0, 1, 2, . . . 
These trajectories undergo 2p+l Andreev reflections and 
2p reflections at tunnel barriers They accumulate 

an action phase 

S yJ = p(-n + + Et tt i) + 2Et lsl - (tt/2 - <p/2). (23) 

One should substitute <p ~ * ~ 4> when interchanging seg- 
ments si and s3, though the relative sign between tt and 
4> does not affect the final result. Here, ttj gives the 
duration of the loop [the sequence between bracket in 
Eq. (f2"2"|) ] and t lsl is the duration of the segment 7 Sl as 
shown on Fig. [5] Class II trajectories split into two sub- 
classes, defined by the following two sequences 



(p) 

Ilia 



(p) 

Tub 



In 
p X 
(e) 



(24a) 



7l2 ) +7l 2 l) +7^+7^ ) +7li ) +T< 1 e) 



(h) 
Tsi 



ls2 



(e) 
'7,2 



(24b) 



This expression sums over all pairs of classical trajecto- 
ries 7I and 72 with fixed endpoints yo and y' , which 



They undergo 2p (IIA) and 2p + 2 (lib) Andreev reflec- 
tions, Ap (Ha) and Ap + 2 (lib) reflections at tunnel bar- 
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riers [40| , and accumulate action phases 

S-yjia = p(-ir-(j) + Et ma ) + (E F + E)t^ n ,(25a) 
S 7 jib = p(— n + (j> + Et£ t ii h ) - (E F -E)t Jn 

+2E(t lsl +t Ja2 ) - (tt-0). (25b) 

Here t^n a and t^m, give the duration of the Andreev 
loops [the two sequences between bracket in Eqs. (|24])]. 
We already see that at E = and (f> = ir, the phase dif- 
ference accumulated by any two members (with different 
p) of a given family vanishes, so that all pair of trajecto- 
ries within a given family resonate. There is however no 
resonance between members of different families. 

In normal chaotic billiards, the stability of periodic 
orbits decreases exponentially with the number of times 
the orbit is traveled [Jl[. The situation is fundamentally 
different in presence of superconductivity, where Andreev 
reflections refocus the dynamics. The stability of a tra- 



jectory is then only given by the product of the stabilities 
along the primitive segments, independent of p. We will 
take these primitive segments as 7 s i and 7 S 3 for class 
I, 7 S i, 7„ and 7 s2 for class II, keeping in mind, how- 
ever, that the contributions arising from the trajectories 
of class I and Ha with p = are more stable than those 
with p > 1, because they do not travel on 7^3 (class I) 
nor on 7 s i and 7 s2 (class Ha). Their contribution is thus 
underestimated in our approach, and our final results, 
Eqs. (j3"Tj) and (|32[) . underestimate the conductance by a 
subdominant correction. 

This enhanced stability applies to trajectories whose 
Andreev loop is shorter than twice the Ehrenfest time, i.e. 
the time beyond which an initially narrow wavepacket 
can no longer fit inside a superconducting lead 0, 0j E3] • 
The relative measure of trajectories of class I (Pj ) and II 
(P/j) is thus 



-P/(te/t d .s) = t ds / / d *«3 exp[-(t s i +t s3 )/T DtS ] = 1 - [1 + t b /td,s] exp[-T B /TD,s], (26a) 

Jo JO 

PTE rTE — t 3 i rTE — t 3 i—t n 

-P//(te/ t d,s) = t d | / di si / dt 'i / dt s2 exp[-(t s i +t„ +t s3 )/r Dj s] 

Jo Jo Jo 

= 1 - [1 + te/td, s + (t e /t d , s ) 2 /2] cx P [-te/t Di s] . (26b) 

We are now ready to evaluate the dominant contributions to conductance close to resonance. We focus on the low- 
temperature regime with E = 0. We start from Eq. (|20[) and pair trajectories by class, noting that for a given class, 
all trajectories have the same stability but differ only by the number of Andreev reflections and normal reflections at 
the contacts to the normal leads, as well as by the different action phases they accumulate along their Andreev loop. 
The sum over classes is then represented by a sum over primitive trajectories, and we perform the substitution 



A 7l a; 2 [. . . ] 7l , 72 — > £ <T (1 - r N )<w r| +p '+ c [. . . ] 7>PiP , . (27) 



E 

71,72 7— primitive p,p'—0 



The exponents a = 1, b = and c = 1 for class I, a = 2, 
6 = and c = for class Ha and a = 2, b = 2 and c = 2 
for class lib are determined by the number of Andreev 
and normal reflections in Eqs. (|2"2"1) and Reflection 
phases do not appear because each time an electron is 
reflected at a tunnel barrier, a hole is reflected later. The 
argument in the double sum over p and p' in Eq. (|27j) has 
to be multiplied by the p-p'-dependent Andreev phase. 
The double sum over p and p' is easily resummed, 

oo 

^ _ p n )a(p+p')+ fc pP+p'+c e i(p-p') (0-7r) _ ^gj 



p,p'— 



r§ (i-r N ) fc 



i-2r s (i-r N ) a cos[7r 



r!(i-r N )*«- 



We next evaluate ^ 7 B^ by relating it to classical trans- 
mission probabilities [29(. Because the sum runs over 
primitive trajectories whose stability is given by their 



electronic (or hole) component alone, one considers that 
the cavity is opened to four normal leads. The sum runs 
over all trajectories such as those sketched in Fig. [6] 
For class I, these trajectories are made up of two legs, 
both of them starting at yo on the L lead, and we write 
E -B 7 = S 7 i, 7 3 B ji B j3- The le § 7 1 makes a normal an- 
gle 9q with respect to the normal to a cross-section of L 
and ends up in lead Si, while 73 makes a normal angle 
9' = —6*o and goes to lead S 2 . Both incidence position 
and angle on Si and S 2 have to be integrated over. 



We define P(Y,Y ;t)SyS8St as the classical probabil- 
ity to go from an initial position and momentum angle 
Yo = (j/oi^o) 011 a lead to within (Sy,S9) of Y = (y, 9) 
on another lead, in a time within St of t. The sum over 
all primitive trajectories of class I can then be rewritten 
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as 



E 



si = 



= PF 



7, primitive 



dt 8 i 



TE—t, 



si 



s3 



dt«3 

tt/2 



tt/2 



7r/2 



d^ d^ 



'SI JS2 J-tt/2 

x P(Y sl ,Yo;i s i)P(Y s3l Yo;t S 3) 

x cos6ocos8' 5{6 + 6' )6(y -y' ). (29) 

The factor ppcos^o measures the injected current, and 
since electrons arc reflected (and not injected) at the L 
tunnel barrier, one has an additional cos#q (instead of 
Pf cos#q). For a given cavity, P is a sum of S- functions 
over all possible classical trajectories. Instead we con- 
sider the distribution averaged over a mesoscopic ensem- 
ble of similar cavities or a small energy interval. This 
gives a smooth function 



(P(Y ri ,Y;t)> = 



4(rw N + rw s ) 



exp[-Vr D ,s] , (30) 



where Wn,s = ^Vn,s7t/A;f gives the width of the normal 
and superconducting leads. We insert ([50]) into (|29p. 
and perform the integrals. Combining the result with 
Eqs. (|2"D1) . (|2"T|) and (|28p delivers the resonant semiclassi- 



cal contribution to T^£, 



0T£>r = 



iVs 



V 2r N A^ N + 2r s iVs 

x (l - (1 + T E /r D;S )exp[-T E /T DiS ] 



(31) 



i-2r s (i-r N )cos[7r 



r 2 (i-r N )s 



The semiclassical contribution to Tjj,^ arises from class 
II trajectories. The geometry is detailed in Fig. [6J and a 
similar calculation as above delivers 



7r 2 r 2 N 2 

\ 1 RLh- — 



N s 



(32) 



87V S V2r N iv N + 2r s A^s 
1 - (1 + te/td.s + t|/2td jS ) exp[-T E /r D ,s] 

i + r 2 (i-r N ) 2 



i-2r s (i- r N ) 2 cos[7r- 



r 2 (i-r N )4' 



The semiclassical contribution to the conductance is 
given by the sum of Eqs. (f3"Tj) and ([52]) . It clearly ex- 
hibits the functional dependence of resonant tunneling, 
where the resonance is however always at the Fermi level, 
and is achieved by setting the phase difference between 
the two superconductors at (f> = tt. This resonance con- 
dition is the same for all semiclassical contributions we 
just discussed. This is why the resonance is macroscopic, 
giving a conductance G oc Nn and not of order one, as is 
usually the case for resonant tunneling. We also see that 
the effect disappears if the superconductors are poorly 
connected to the normal cavity, Ts — > 0, as should be. 





Figure 6: The effective four-terminal cavity used in the cal- 
culation of B~ , and the geometry of class I (top panel) and 
class II (bottom panel) trajectories. 



Comparing Eqs. (J31J) and (J32J) to Eq. (JT3J) we see 
that semiclassical contributions matter only close to res- 
onance. In the tunneling regime, they increase the con- 



ductance by a factor 



from G(tt)/Go oc T n N^ 



to G(tt)/Gq oc ATn, resulting in a peak-to-valley ratio of 



the conductance G((j> = ir)/G(4> = 0) oc r r 



In most 



instances, (T£^(7r)) r » {T^ c L (n)) r , and the sharpness of 
the resonance peak, measured by its width at half height, 
is proportional to F. 



IV. NUMERICAL SIMULATIONS 



We compare our theory with numerical simulations on 
a quantum mechanical model of the Andrccv billiard, the 
Andrccv kicked rotator 11211. 
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The Andreev kicked rotator 



with 



One-dimensional maps give a stroboscopic description 
of dynamical systems. An example is the kicked rotator, 
describing the dynamics of a particle confined to a circle 
which is kicked periodically in time with period To and a 
kicking strength depending on its position Q . The quan- 
tized version of this map is characterized by its Floquet 
operator, the unitary operator giving the time-evolution 
from one kick to the next, i.e. from time t to t + To- We 
represent the Floquet operator as a M x M matrix with 
matrix elements (from now on we measure times in units 

Of T = 1) 



F„ 



-(z7r/2M)(n 2 +m 2 ) 



(UQU^) r 



jVf -1 / 2 e (2iri/M)nm 



-(iMK/2-ir) cos(27rn/M) 



(33a) 
(33b) 
(33c) 



Here K is the kicking strength. For K > 7 the clas- 
sical dynamics is chaotic, with classical Lyapunov ex- 
ponent A ~ In [A/ 2]. The eigenvalues of F are e l£m 
(to = 1...M), and the quasi-energies e m of the closed 
system (without any leads) have an average level spacing 
5 = 2n/AI. The dynamics of the electron is determined 
by P, while the time evolution of holes is given by F*. 

The kicked rotator model describes a closed system. 
We open it by inserting leads in a subspace of the Hilbcrt 
space, parametrized by the indices {mi, TO2, . . . TO2at n } 
for the normal leads and {l±, I2, ■ ■ ■ l2N s } f° r the super- 
conducting leads. The coupling to the leads is described 
by a 2 An x M (2N s x M) projection matrix P N (P s ) 
with matrix elements 



[pNl]ram 
[Ps]nm 



sJT^ ifm = n G {mi, . . .to 2 at n }, 
otherwise, 

v/Ts if m = n G {h, ■ . ■ hN s }, 
otherwise. 



(34a) 
(34b) 



The average time between Andreev reflections equals 
t d.s = -^^/2rgAg. In our numerical investigations we 
restrict ourselves to Ts = 1. The Ehrcnfest time of the 
closed Andreev billiard reads r E = A" 1 In (Af/M) [gL fl^| . 

We perform two sets of simulations, the first is devoted 
to the spectrum of a closed Andreev billiard, the second 
to the conductance through a chaotic Andreev interfer- 
ometer. 



1. The closed Andreev kicked rotator 

Both electrons and holes have an energy near the Fermi 
energy, and therefore we assume that Andreev reflection 
is perfectly retroreflccting. The Floquet matrix of the 
Andreev billiard is then 



F 
F* 



(35) 



3>s = 



1 - PgPs 



TJ$> 



-iPge 



Ps 



1 



PlPs 



(36) 



The coupling between electrons and holes is modelled 
by ^Ps, which induces Andreev reflection for electron or 
holes touching the superconducting leads, i.e. induces 
transitions between the components {l±, . . . hN s } °f the 
wavefunction in the two quasiparticle sectors. The phase 
difference between the two superconductors is contained 
in the 2N$ x 2N$ matrix with elements S nm (j)/2, for 
n < As and —S nm 4>/2 for n > Ng. The matrix 5" can 
be diagonalized efficiently using the Lanczos technique in 
combination with the Fast-Fourier- Transform algorithm 
[42I ]. This allows to reach the large system sizes needed 
to explore the regime of finite Ehrenfest time. 



2. The kicked Andreev interferometer 

We first write the (2A N + 2A S ) x (2A N + 2A S ) energy- 
dependent scattering matrix of our four-terminal struc- 
ture (without superconductor). It is related to the Flo- 
quet operator via [43| 



5(e) = -(1 -PP 



T\l/2 



+ P- 



1 



F(l - P T P) 



FP 1 



The (2A S + 2 An) x M coupling matrix 



P 



Pn 
P s 



(37) 



(38) 



describes coupling to all four leads. The four-terminal 
scattering matrix S is made of sixteen sub-blocks Sij, 
where the subindices i,j = L, R, SI and S2 la- 
bel either normal or superconducting terminals. From 
these sub-blocks, the exact expression for the conduc- 
tance, Eq. (fl6)) . can be computed, following Ref. [25| . 
The transmission amplitudes giving the probabilities in 
Eq. p7|) can be written as 



.ho 



fl-tfic'fi^l + cnc'O*)^, (39a) 
:-ifc*fi*(l +cClc*n*)- 1 d, (39b) 



in terms of the matrices 

Sl,l Sl.r 
Sr,l Sr : r 

Ssi,Sl Ssi,S2 
SS2,S1 5*52,52 

e"* 2 
e^/ 2 



b = 



Sr,si 

Ssi,L 
Ss2,L 



Sl,S2 
Sr,S2 

Ssi,R 
SS2,R 



n 



For numerical stability and efficiency, the transmis- 
sion matrices r of Eqs. (|39[) are constructed from self- 
consistent linear equations for the output amplitudes as 
a function of unity input on each An channel, one at a 
time, and not via direct inversion. 
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1 2 3 1 2 

E/Ej s E/5 

Figure 7: Left panel: density of states for the Andreev kicked 
rotator with (ft = 0, M = 2048, iVs = 25, K = 147, r s = 1. 
The solid line gives the prediction of Eqs. Q and Q. Data 
are averaged over 50 different lead positions. Right panel: 
density of states for the Andreev kicked rotator with (ft — tt, 
M = 2048, iV s = 25, T s = 1. The solid line gives the predic- 
tion of Eq. ((SJ). Data are averaged over 50 different lead posi- 
tions and 50 values of the kicking strength K £ [144.5, 147.5] 



B. Numerics in the universal regime 

We start with a comparison of the analytical predic- 
tions of RMT with numerical results, by choosing the 
parameters of the Andreev kicked rotator in such a way 
that te <C 71) s- In Fig.[7]we show the density of states for 
(f> = and = tt. We compare our results with the the- 
oretical predictions of Eqs. and ([4]) for = 0, while 
for (ft — tt wc use Eq. ((SJ). In both cases we find perfect 
agreement between numerics and RMT predictions. 

We next turn our attention to the conductance. We 
choose an Andreev billiard with the same parameters as 
in Fig. but couple it to normal leads. The conductance 
as a function of the phase difference is shown in Fig. [5] 
The numerical results, for several values of the param- 
eters Tn and JVn (characterizing the normal leads), are 
compared with the circuit theory prediction of Eqs. pip 
and (fT5|) . The agreement is excellent. In the tunnel- 
ing limit Tn — > the conductance is strongly suppressed 
for 7^ tt while it is restored to its value in absence of 
superconductor when = tt, as predicted by Eq. (|13p . 

In Fig. [21 we show the conductance clS cX function of 
Tn, both for (ft = and (ft = tt. For = the circuit 
theory of Eqs. (fTTj) and (fT2"|) (solid lines) predict a curve 
that depends only on the ratio a = Ns/Nr- The nu- 
merical data line up on the analytical curves. For (ft = tt 
the analytical prediction of circuit theory coincides with 
the conductance of the billiard in the absence of super- 
conductors i.e. G/Go = Tn-/Vn/2. As T^N^ is reduced, 
however, the data deviate from this prediction. We at- 
tribute this to the fact that the circuit theory results are 
accurate to leading order in T-^N^ only. This argument 
is corroborated by the fact that the data for the largest 
A r N remain close to the circuit theory prediction down to 
lower values of Tn • 

We next present in Fig. [TO] data for the conductance 




Figure 8: Conductance of the Andreev kicked rotator as a 
function of the phase difference (ft between the two supercon- 
ductors. For all curves M = 2048, iV s = 25, K = 147. The 
closed symbols correspond to TVn = 100 and the open sym- 
bols to TVn = 50, with Tn = 1 (squares), Tn = 0.2 (circles), 
Tn = 0.1 (triangles) and Tn = 0.01 (diamonds). The solid 
lines are the relevant predictions from circuit theory. Every 
data point is an average over 200 realizations. 




Figure 9: Conductance of the Andreev kicked rotator as 
a function of T N for K = 147, M = 2048, r s = 1, at 
(ft = (closed symbols) and (ft — tt (open symbols), and 
N N /N S = 100/50 (squares), 100/25 (circles), 100/12 (upward 
pointing triangles), 200/25 (diamonds) and 400/50 (down- 
ward pointing triangles). Every data point is an average over 
50 realizations. The solid lines are the predictions from circuit 
theory for (ft = and a = N s /N^ = 0.5,0.25, and a = 0.12, 
while the dashed line corresponds to (ft — tt. 



at finite excitation energy. For (ft = tt, the conductance 
stays at its universal, E = value, G = FnA^n/2, in- 
dependent of E. For = however, it increases from 
G = T^N^{l + 2N^/N s ) to G = r N N N /2 as E increases. 
The crossover occurs on the scale ~ 0.6i?T,s of the = 
Andreev gap (see the top panel of Fig. [7]), which illus- 
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Figure 10: Transmission as a function of the normalized ex- 
citation energy E/E T ,s for M = 2048, N N = 100, K = 147, 
Tn = 0.05, at <j> = (black symbols) and cj> = tt (red symbols) 
for Ns = 25 (triangles), 50 (diamonds), 100 (squares) and 200 
(circles). The black curve shows the universal conductance at 
d> = tt for E = 0. 



to 
w 
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Figure 11: Density of states of closed Andreev billiards, for 
<j> = (left panel) and <j> = tt (right panel), in the universal 
regime with M = 2048, N s = 25 and K = 148 (red curves, 
Te/td,s = 0, data averaged over 2500 samples) and in the 
semiclassical regime with M = 8192, Ns = 800 and K = 10 
(green curves, te/td,s = 0.55, 500 samples) . The black 
curves indicate the universal predictions [J, [35|]. Note the 
different horizontal and vertical scales. 



trates further the connection between transport through 
the Andreev interferometer and the spectrum of the cor- 
responding closed Andreev billiard. From the data in 
Fig. [TU1 we do not expect a fundamental departure from 
the zero-temperature theory in the universal regime, as 
long as the temperature lies below the Thouless energy, 
T < £t,s- The situation is different in the semiclassical 
regime, which can already be inferred from the Ns = 200 
data presented in Fig. [TU1 Despite a rather small ratio 
Te/td,s ~ 0.13, this set of data already exhibits a peak 
at E = above the universal conductance for 4> = tt. 
This peak however quickly disappears with E <C -E»r,S- 

We finally stress that for the data presented in Figs. [7] 
HI and [HI t e /t d ,s < 0.02. We found that, to get a good 
agreement between our numerics and RMT predictions 
at <f> = 7T, it is necessary to set te/td.s to much smaller 
values than in our earlier numerical works with = [l2l . 
0, [3 ■ This is especially true for the density of states at 
4> = tt (where we needed to go down to te/td,s ~ 10 -3 ), 
and is already a good indication that the semiclassical 
effect at </> = tt is much larger than the manifestations of 
Ehrcnfest physics investigated so far. 



C. Semiclassical limit 

We next investigate the spectrum of closed Andreev 
billiard at finite te/td,s- Density of states of closed bil- 
liards are shown in Fig. [TTJ The density of states is only 
weakly dependent on te/td,s at <f> = (left panel). At 
4> = tt however, the semiclassical regime sees the emer- 
gence of a giant peak in the density of states at the Fermi 
energy (right panel). 

In Fig. 1121 we next show a conductance resonance curve 
in the semiclassical regime. We obtain very good agree- 



ment between the numerical data (circles) and the ana- 
lytical prediction (green solid line) with te/td,s — 0.79. 
Without the semiclassical contribution, the universal pre- 
diction of Eqs. (|lH412p in the full range € [0,2tt] is 
given by the red line. While the latter fits the resonance 
curve far away from resonance, the agreement breaks 
down close to (f> = tt, where universal contributions give 
a prediction G(tt) = 2Gq, too small by an order of mag- 
nitude. The left inset in Fig. [TJ] illustrates the increase 
of the peak height and narrowness as the semiclassical 
parameter M increases. The four sets of data in this in- 
set correspond to a fixed set of classical parameters, with 
the electronic wavelength decreasing by factors of four 
from one curve to the next, starting from the bottommost 
(blue) curve. The conductance increases at each step be- 
cause the number of conduction channels scales linearly 
with M. In absence of semiclassical contributions, these 
four curves would exhibit the same peak-to-valley ratio, 
but here they do not. This is quantified in the right inset 
to Fig. [TJJ where we show both the peak height and the 
peak-to-valley ratio corresponding to the same classical 
configuration as in the main plot, while varying M. 

We can connect the data in Figs. [TT1 and [T2l since they 
correspond to the same closed Andreev billiard. The 
broadening of levels due to the external normal leads is 
SE ~ rN-/VN<5/27r. For the resonance in Fig. [TJl this 
gives 8E = 26 /tt. From the green curve in Fig. 111! we 
estimate that this value of SE covers about 20 to 30 lev- 
els around the Fermi energy. Those states condense on 
classical orbits bouncing on average M/2N$ = 5 times 
at the billiard's boundary. At each bounce, they have a 
probability ~ 1 — 2N-^/(M — 2N$) not to touch a nor- 
mal lead. Therefore, since [1 - 2N N /(M - 2N S )} 5 « 1/2 
for the data in Fig. [TJ] with A^n = 400 channels, we es- 
timate that about half of these 20-30 states couple to 
the leads. One thus expects a resonant contribution to 
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Figure 12: Conductance through a chaotic Andreev in- 
terferometer vs. the phase difference <j> between the two 
superconductors. Circles are numerical results obtained from 
the Andreev kicked rotator with T N = 0.01, M = 8192, 
M/JV N = 20, M/N s = 10 and K = 10. The green curve is the 
analytical prediction obtained by summing the semiclassical 
resonant contributions with the universal prediction. The 
red curve gives the universal prediction obtained from circuit 
theory. Left inset: Numerical data for the same classical 
parameters M/N N - 20, M/Ns = 10 and K = 10 as in the 
main plot, for M — 128 (blue curve), 512 (green curve), 2048 
(red curve), to 8192 (black curve). Note the change in peak- 
to-valley ratio. Right inset: peak-to-valley ratio G(n)/G(0) 
(black squares) and peak conductance G(n) (red circles) as a 
function of M controlling the quantum-to-classical crossover, 
for the same classical configuration as in the main plot. Data 
are averaged over 150-1000 sample realizations. 

the conductance ~ 10 — 15Go, which agrees quite well 
with the numerical value of G(ir) ~ 16Go for the data of 
Fig. [T2l with a universal contribution of 2Go- 

In Fig. [13] we show the conductance as a function of the 
number of channels in each superconducting lead Ns , for 
two different values of Tn- We see that the conductance 
depends strongly on TVs, in contrast to the RMT predic- 
tion G = GoTnA^n/2 which is independent of it. In the 
parameter range of the figures, the contribution of Eq. 
(|32| can be neglected and the contribution of Eq. (j3l"j) 
dominates. For <p = n and Tn — > 0, it can be simplified 

to T'H = (7tiV N /16)(l-(l + T E /r D ,s)exp[-T E /rD,s])- 
This contribution depends strongly on N§ since both the 
Ehrcnfcst time and the time td,s between Andreev re- 
flections depend on N$. It is therefore the relative mea- 
sure of trajectories of class I which causes the order-of- 
magnitude increase of conductance with N$ observed in 
Fig. [13] The solid lines give the theoretical prediction 
obtained by summing the semiclassical resonant contri- 
bution, Eqs. (|3ip and (|32p , and the universal contribution 
obtained from Eqs. (TlTj) and (TT2"|) . There is good agree- 
ment between the numerical data and the theory as long 
as the superconduting dwell time is long enough. For 
Ng > 400, td,s < 5, and we suspect that short processes 
across the billiard increase the conductance well above 



Figure 13: The conductance of the Andreev kicked rotator as 
a function of N s for K = 10, M = 4096, r s = 1, <j> = tt, An = 
200 and Fn = 0.03 (closed circles), Fn = 0.01 (open circles). 
Every data point is an average over at least 80 realizations. 
The lines represent a sum of the semiclassical prediction Eq. 
([3Tj) in the tunneling limit and the RMT result of Eq. JT3J for 
Tn = 0.03 (solid red line) and Fn = 0.01 (dashed green line). 



the analytical prediction. 



o 




1 ^ ■ ■ ■ ■ ^ 1 

100 200 300 400 500 600 700 
N N 



Figure 14: The conductance of the Andreev kicked rotator 
as a function of A N for K = 10, M = 4096, A s = 400, 
Fs — 1, 4> ~ tt and Fn = 6/An (closed circles), Tn = 2/An 
(open circles). Every data point is an average over at least 
80 realizations. The lines represent a sum of the semiclassical 
prediction Eq. f|31f) in the tunneling limit and the universal 
RMT result of Eq. JTSJl or T N = 6/A N (solid red line) and 
Tn = 2/An (dashed green line). 

Finally, we show on Fig. [T2] the conductance as a func- 
tion of JVn. The product T^N^ is kept constant, so that 
any observed N-^ dependence arises from nonuniversal ef- 
fects. Even though we always have moderate Ehrcnfcst 
times te < td.s, the numerical results are very different 
from the constant conductance predicted in the universal 
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regime. The semiclassical contributions dominate and we 
see an increase of G with 7V N . Here, te/td,s — 0-6 is con- 
stant, and the increase of the conductance is predicted 
to be linear, see Eqs. (f3"Tj) and (j3"2")h There is reasonable 
agreement between theory and numerical data, with sig- 
nificant deviations in the regime of small TVn, where we 
leave the tunneling regime. 

We conclude that our numerical data fully confirm our 
analytical results. 

V. CONCLUSION 

We have presented analytical and numerical investiga- 
tions of spectroscopy and transport in quantum chaotic 
Andrccv billiards coupled to two superconductors. We 
identified two regimes with very different behaviors. In 
the universal regime, the presence of the superconductor 
opens up a gap in the density of states of the Andreev bil- 
liard at the Fermi level, which reduces the conductance 
through the system once it is attached to two external 
normal leads with tunnel contacts. This gap closes when 
the phase difference between the two superconductors is 
(p = 7T, and in the tunneling limit, the normal conduc- 
tance that the system would have in the absence of su- 
perconductivity is recovered. As the Fermi wavelength 
is reduced, Xf/Ws, Xf/L c <C 1, one enters a new, semi- 
classical regime, where the situation is only marginally 
different unless the phase difference between the two su- 
perconductors is close to 7r. When this is the case, we 
found an order-of-magnitudc enhancement of the tunnel- 
ing conductance, and identified the mechanism behind 
this enhancement as resonant tunneling through a macro- 
scopic number of quasi-degenerate levels at the Fermi en- 
ergy of the corresponding closed Andreev billiard. 

Our calculation focused on the zero-temperature limit 
of transport, and we finally comment on the effect of a 
finite temperature. From Fig. 1101 we concluded that, 
in the universal regime, the shape of the conductance 
oscillations do not differ from those presented in Fig.[8l as 
long as the temperature lies below the superconducting 



Thouless energy, Et,s = TsNsS/2ir. The semiclassical 
contribution to the conductance near <f> = tt is however 
much more sensitive to changes in temperatures than the 
universal one. Finite excitation energies indeed lead to 
an additional action phase accumulated along transport 
orbits. The expressions in Eqs. (J^U) and ([25)1 are now 
combined into the following expression 

J (J _ p n )a(p+p')+b pp+p'+c e i(p-p')(0-' r ) (40) 

r D,s Pi j/ = o 

x T E di sl I™'** 1 dt s3 e-^ +t ^ T ».s +iE( *- p 'V H \ 
Jo Jo 

The accumulation of an E'-depcndent action 
phase results in a cut-off for the summation over 
p and p' , effectively suppressing contribution with 
p — p' > PmaL = (h/E td,s) 2 - One ends up with a 
prefactor 1 — (1 — rN) p ™ a>x multiplying Eqs. (|3Tj) and 
(|32[) . which suppresses the semiclassical contribution 
already for temperatures significantly smaller than 
£/r,s- At finite temperature, the resonant conductance 
enhancement discussed above will thus be maximal for 
an intermediate value of Tn optimizing macroscopic 
resonant tunneling while minimizing thermal averaging. 
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